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Abstract 

Accounting for the two independent correlation functions of the QCD vac- 
uum, we improve the simple and consistent description given by the model of 
the stochastic vacuum to the high-energy pp and pp data, with a new determi- 
nation of parameters of non-perturbative QCD. The increase of the hadronic 
radii with the energy accounts for the energy dependence of the observables. 

PACS Numbers : 12.38 Lg , 13.85 -t, 13.85 Dz , 13.85 Lg . 
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1. Profile Function for Hadron-Hadron Scattering 

In the effort to explain the properties of soft high-energy scattering in a non-pertubative 
QCD framework, we reconsider the applicatioJ of the model of the stochastic vacuum 
of non-perturbative QCD to hadron-hadron scattering. This approach combines the pa- 
rameters describing properties of the QCD field (gluon condensate, correlation length) with 
those describing the colourless hadrons. The model presents the characteristic features of the 
pomeron exchange mechanism of Regge phenomenology, where vacuum quantum numbers 
are exchanged between hadronic structures. For all hadronic systems the total cross-sections 
increase with the energylli somewhat between s - 0808 and log 2 (s/so) within the present ex- 
perimental range, and this is explained in terms of the energy dependence of the effective 
hadronic radii. 

Several models relate the total high-energy cross-sections and slope parameters to the 
hadronic radii@i. This is a characteristic feature also of the model of the stochastic vac- 
uum which gives specific predictions for the size dependence of the observables for different 
hadronic systems^ The knowledge of the hadronic structures required for the description 
of the soft high-energy data does not go beyond the information on their sizes, the sim- 
plest transverse wave-function giving all information required for the determination of the 
observables. 

The model of the stochastic vacuum deals successfully with non-perturbative effects 
both in low-energy hadron physics! and in soft high-energy scattering^ . The treatment of 
scattering is based on the concept of loop-loop interaction, which allows a gauge-independent 
formulation for the amplitudes. The loops, formed by the quark and antiquark light-like 
paths in a moving hadron, have their contributions added incoherently, with their sizes 
weighed by transverse hadronic wave-functions. 

The principles and methods used for the evaluation of the observables of high-energy 
scattering in the model of the stochastic vacuum have been fully described before^. In the 
present work we present the results of a more complete calculation in which both Abelian 
and non- Abelian contributions to the field correlator are taken into account. The dynamics 



is generated by correlations in the QCD vacuum field, according to the pioneering work 
of NachtmannS. With the assumption that the correlator is independent of the reference 
point of the gauge field, its most general formi contains two tensor structures, each with a 
correlation function, with a weigth parameter k between them. 

To write the expression for the profile function for hadron-hadron scattering, we introduce 
the notation R(I, J), where the first index (1=1,2) specifies the loop, and the second specifies 
the particular quark or antiquark (J=l or 2) in that loop. Fig.l shows a projection of the 
collision of two mesons on the transverse scattering plane. The vectors Q{K, L) in the 
transverse plane connect the reference point C to the positions of the quarks and antiquarks 
of the loops 1 and 2. The quantity ip{K, L) is the angle between Q(l, K) and Q(2, L). 

The eikonal function of the loop-loop amplitude is written 

X(b, R(l, 1), R(2, 1)) = k[ - cos V>(1, 1) J[Q(1, 1), Q(2, 1), 1)] 

-cos^(2,2) J[Q(1,2),Q(2,2),^(2,2)] 
+ cos^(l,2) J[Q(1,1),Q(2,2), 0(1,2)] 
+ cosV(2,l)/[Q(l J 2),Q(2,l) > V(2,l)] 
+(!-«)[- W[Q(1, 1), Q(2, 1), ^(1, 1)] - W[Q(1, 2), Q(2, 2), -0(2, 2)] 

+W[Q(l,l),g(2,2),^(l,2)] + ^[g(l,2),g(2,l),^(2,l)]] . (1) 

The quantities / which represent the non-Abelian contributions are given by integrations 
along the dashed lines of the figure 

I[Q(l,K),Q(2,L)MK,L)] = ^[^j 

x{Q(l,K) / [Q(l,K) 2 + x 2 -2xQ{l, K) cos i;{K,L)} 
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Q(2, L) 2 + x 2 - 2xQ{2, L) cosip{K, L) 



dx} 
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with Q(K,L) = \Q(K,L)\. The quantities W, which come from the non-confining part of 
the correlator, are given by 

W[Q(l,K),Q(2,L)^(K,L)} = g 2 ^ 



x [Q(l, K) 2 + Q(2, L) 2 - 2Q(1, L) cos V>(#, L)] 



3/2 



^7g(l, if) 2 + Q(2, L)2 - 2Q(1, tf)Q(2, L) cos^(K, L) 



K 3 

From the eikonal function x we form the loop-loop amplitude JL 1 L 2 (b, Ri, R2) through 



(3) 



J LL ,(b, R x , R 2 ) = *J L l>& R x , R 2 ) = (4) 

v ' ' ; [(g 2 FF)} 2 144-576 v; 

where Ri and R 2 are shorthand notations for R(l, 1) and R(2, 1) respectively. 

The hadron-hadron amplitude is constructed from the loop-loop amplitude using a simple 
quark model for the hadrons, by smearing over the values of R x and R 2 with transverse 
wave-functions ip(R). Taking into account the results of the previous analysis of different 
hadronic systenJ, in the present calculation we only consider for the proton a diquark 
structure , where the proton is described as a meson, in which the diquark replaces the 
antiquark. Thus these expressions apply equally well to meson-meson, meson-baryon and 
baryon-baryon scattering. 

For the hadron transverse wave-function we make the simple ansatz 



MR) = pl^ir ex P (~ R l s n) > ( 5 ) 

where Sh is a parameter for the hadron size. We then write the reduced profile function of 
the eikonal amplitude 

J HlH2 (b,s h s 2 ) = J rf 2 i?i J d 2 R 2 J LlL2 (b,RuR2) \MRi)\ 2 \MR2)\ 2 , (6) 

which is a dimensionless quantity. The dimensionless scattering amplitude is given by 

T HlH2 = is{(g 2 FF)a 4 ] 2 a 2 J d 2 b exp {iq ■ b) J HlH2 (b, S x , S 2 ) , (7) 

where the impact parameter vector b and the hadron sizes S\ , S 2 are in units of the 
correlation length a, and q is the momentum transfer projected on the transverse plane, in 



units of 1/a, so that the momentum transfer squared is t = — \q\ 2 /a 2 . For convenience we 

have explicitly factorized the dimensionless combination (g 2 FF)a 4 . The normalization for 

Th 1 h 2 is such that the total and differential cross-sections are given by 

1 da el 1 

a T = - Im T HlH2 , — — = — — \T HlH2 \ 2 . (8) 
s at lbirs^ 

For short, from now on we write Jib) or J(b/a) to represent J(b, Si, S 2 ). To write 

convenient expressions for the observables, we define the dimensionless moments of the 

profile function (as before, with b in units of the correlation length a) 



d 2 bb k J{b) , = 0, 1,2, ... (9) 

which depend only on Si/ a, S^/a, and the Fourier-Bessel transform 

I{t) = J d 2 b J {basJ\t\) J{b) , (10) 

where Jo(ba^J~\t~\) is the zeroth-order Bessel function. Then 

T HlH2 = iS [(g 2 FF)a 4 ] 2 a 2 I(t) . 

Since Jib) is real, the total cross section a T and the slope parameter B (slope at t = 0) are 
given by 



' 'V .Ka>. (11) 



2 „2 _ r^„2 

t=o 2 J 



The model conveniently factorizes the QCD strength and length scale in the expressions 
for the observables, and the correlation length appears as the natural unit of length for 
the geometric aspects of the interaction. These aspects are contained in the quantities 
Io(Si/a, Sz/a) and ^{Si/a, £2/0), which depend on the hadronic structures and on the 
shapes and relative weights (parameter k) of the two correlation functions. It has been 
shown beforeB that for the case « = 1 (pure confining correlator) the two moments have 
simple form as functions of S/a, and in the present work we have obtained similarly simple 
expressions for any value of k. We concentrate on the range about k = 3/4, which is 
indicated by lattice results. 



It is important that the high-energy observables a T and B require only the two low 
moments Jo , 1% of the profile functions. We observe that in the lowest order of the correlator 
expansion the slope parameter B does not depend on the value of the gluon condensate 
(g 2 FF) and, once the proton radius S is known, it may give a direct determination of the 
correlation length. 

2. pp and pp systems and QCD Parameters 

In order to have a wide range of data to extract reliable information on QCD parameters, 
we concentrate here on the pp and pp systems, for which Si — S% — S. The curves for 
I = a T I [{g 2 FF) 2 a 10 } and K = B/a 2 can be parametrized as functions of S/a with simple 
powers, with good accuracy. The convenient expressions are 



The numerical values for the parameters are given in table 1. 

The parametrizations for a T and B are very convenient for comparison of the results of 
the model of the stochastic vacuum with experiment. In the present analysis we take into 
account all available data on a T and B in pp and pp scattering , which consist mainly!'! 
of ISR (CERN) measurements at energies ranging from y/s = 23 GeV to \fs = 63 GeV, 
of the a/s = 541 — 546 GeV measurements in CERN SPS and in Fermilab, and of the 
\/s = 1800 GeV information from the E-710 Fermilab experiment. Besides these, there is a 
measurement^ of o T = 65.3 ± 2.3 mb at ^/s = 900 GeV and there are the measurements of 
a T = 80.6 ± 2.3 and B = 17.0 ± 0.25 GeV~ 2 in Fermilab CDF0 at s/s = 1800 GeV which 
seem discrepant with the E-710 experiment at the same energy. A measurement by Burq et 
al@ at y/s = 19 GeV, which was taken as input in a previous calculation^, seems to disagree 
with the ISR data, presenting a too high slope B = 12.47 ± 0.10 GeV -2 . 

In fig. 2 we plot o T and B against each other. At the ISR energies we use <rJ om = 
(21.70 mb)s 0808 and the pp slope parameter as representative of the non-perturbative con- 
tributions which are the concern of the model of the stochastic vacuum. At the highest ener- 




(12) 
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gies (541-1800 GeV) the process is essentially non-perturbative and dominated by pomeron 
exchange. The relation between the two observables is fitted with the form 

B = B A + C A (a T ) A . (13) 

This form, which is represented by the solid line in the figure, is suggested by the results of 
the calculations with the model of the stochastic vacuu where an interpretation for the 
meaning of the parameters is given in terms of QCD and hadronic quantities. 

The parameters in the model that are fundamentally related to QCD are the weigth 
k, the gluon condensate (g 2 FF) and the correlation length a. The hadronic extension 
parameter Sh accounts for the energy dependence of the observables. We here show how 
these quantities can be evaluated using only high-energy scattering data. 

The proton radius can be eliminated from eqs(|ll]) and fll2"D , and we obtain a relation 
between the observables a T and B at a given energy 

{B ~ Va ) = (< g 2 FF > a 4 ) 2 *//* {~^J ' (M) 
This form is the same as given by eq.(^), with an obvious correspondence of parameters. 
We first remark that the exponent A = 8/(3 does not depend on QCD quantities and is 
almost constant (equal to about 0.75) in the range of values of k that are obtained in lattice 
calculations (k ~ 3/4). Then we fix A = 8/(3 = 0.75 and are left with only two free 
quantities in both energy independent relations (|T3|)and (|T4j ), and they can be determined 
using as input the two clean experimental points for a T and B at 541 and 1800 GeV. The 
results for the parameters are nearly the same as obtained in a free fitting of all points. 
With a T in mb we have 

B A = r]a 2 = 5.38 GeV" 2 = 0.210 fm 2 , C A = 0.458 GeV" 2 . (15) 

Of course these results are subject to uncertainties. We have adopted an ansatz for 
the correlation function, which is arbitrary (although numerically it could not have very 
different shape). There is some uncertainty also in the determination of the parameters 



a,/3 . . . representing the final results of the numerical calculation. On the other hand, the 
model gives a rather unique prediction for A = 8/(3 = 0.75 and this value is well sustained 
by the data as shown in fig. 2 . 

We have made calculations with several values of k, but to be specific, we borrow from 
lattice calculation the value k = 3/4, and then use as parameter values the numbers shown 
in table 1. Taking into account the experimental error bars in the input data at 541 and 
1800 GeV, we obtain 

k = |, a = 0.32 ± 0.01 fm , < g 2 FF > a 4 = 18.7 ± 0.4 , < g 2 FF >= 2.7 ± 0.1 GeV 4 . 

(16) 

With the value k = 33/40 obtained in more recent lattice results^ the central values change 
only slightly to 

k = 33/40 , a = 0.33 fm , < g 2 FF > a 4 = 19.2 , < g 2 FF >= 2.6 GeV 4 . (17) 

The results of the pure SU(3) lattice gauge calculation^! for the correlator 
(FC(x, 0) F£(0, 0)) a have been fitted^ with the same correlation function used in the present 
work. The correlation between the values of < g 2 FF > and a that was then obtained can 
be well represented by the empirical expressions 



1.1122 9 0.01813 oa n~ 

A l = ^TMo- > <£ < FF> =^nm- » < 9 p F > a 4 = 0.0172^ , (18) 



with the lattice parameter in MeV, a in fm, and < g 2 FF > in GeV 4 . This correlation 
is displayed in fig. 3, for values of Al in the usual range. The point representing our results 
in eq.flTB]) is marked in the same figure. The dashed line represents the relation with the 
string tension obtained in the application of the model of the stochastic vacuum! to hadron 
spectroscopy; for our form of correlator, this relation isS 

KifFF) = . (19) 

As can we may see in the figure, the constraints from these three independent sources of 
information are simultaneously satisfied, building a consistent picture of soft high-energy pp 

8 



and pp scattering. The (pure gauge) gluon condensate is well compatible with the expected 
value. The lattice parameter A^ and the string tension p are also in their acceptable ranges. 
As we describe below, the resulting proton size parameter S p takes values quite close to the 
electromagnetic radius0. 

In this model the increase of the observables with the energy is due to a slow energy 
dependence of the hadronic radii. An explicit relation is obtained if we bring into eqs.fllTD 
and fljjD a parametrization for the energy dependence of the total cross-sections, such as 
the Donnachie- Landshoff! form. In this case we obtain for the proton radius 
S (s) - a 1 1 f 21-7mb \ VP //3 

The energy dependence, given by a power 0.0808//3 of s is very slow, and the values obtained 
for S p are in the region of the proton electromagnetic radiusEl, which is R p = 0.862 ± 0.012 
fm. However the Donnachie-Landshoff parametrization for the total cross-sections is very 
convenient to identify the pomeron contributions at the ISR energies, but does not give the 
best representation of the cross-sections at higher energiesi. Using eqs. (|TTD and (|T2| ) and 
directly the data at 541 and 1800 GeV, we obtain the values for the proton radius that are 
shown by the small squares in fig. 4, where a log e scale is used for \/s. It is remarkable that 
we have an almost linear dependence, which can be represented by 

S p (s) = 0.671 + 0.057 log (fm) , (21) 

with \fs in GeV. With this form for the radius, which is shown in dashed line in fig. 4, the 
cross-sections evaluated at very high energies rise with a term log^ y/s, and are smaller than 
predicted by the power dependence of Donnachie-Landshoff. However, since j3 ~ 2.8, they 
still violate the bound log 2 \fs. This may be corrected using a power 2/(3 instead of 1 in the 
parametrization for S p (s), and we then obtain 

S p (s) = 0.572 + 0.123 [log yfsf 72 (fm) . (22) 

This form is shown in solid line in fig. (4). Clearly it gives a good representation for the 
existing data. At 14 TeV, which is the expected energy in the future LHC experiments, 
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we obtain S p (14 TeV) = 1.19 fm = 1.38 R p = 3.7a and the model predictions for the 
observables are a T =92 mb and B=19.6 GeV -2 . The dashed line representing eq. (|2T|) leads 
at the same LHC energy to a T =97 mb and B=20.1 GeV -2 , while the Donnachie-Landshoff 
formula leads to a higher value <t t =101.5 mb. A recent fit of all data! gives a power 
2.25 ± 0.35 in the logarithm, and predicts a T = 112 ± 13 mb at 14 GeV. 

The non-perturbative QCD contributions to soft high-energy scattering are expected to 
be dominant in the forward direction, thus determining the total cross-section (through the 
optical theorem) and the forward slope parameter. The model, as it is presented in this 
paper, leads to a negative curvature for the slope B(t), which decreases as \t\ increases. The 
data however shows an almost zero curvature of the peak, so that above some value of \t\ the 
model leads to too high values of the differential cross-section. Rather small changes in the 
form of the profile function J(b) that enters in the expression for the scattering amplitude 
may modify this behaviour. These changes can be made phenomenologically, introducing a 
form factor^. However, in the present work we wish to keep the fundamental characteristics 
of the model, which is that of a pure QCD based calculation, with a unique set of quantities 
governing all systems at all energies. 

3. Conclusions 

Extending the previous calculation of high-energy observables, now including more data 
in the analysis, and taking into account the two independent correlators of the QCD gluon 
field, we give a unified and consistent description of all data on total cross-section and slope 
parameter for the pp and pp systems, from y/s ~20 to 1800 GeV. The non-perturbative QCD 
parameters determining the observables, the gluon condensate and the correlation length of 
the vacuum field fluctuations, are determined. The third quantity entering the calculations 
is the transverse hadron size, which has a magnitude close to the electromagnetic radius, 
and whose slow variation accounts for the energy dependence of the observables. 

The model allows a very convenient factorization between the QCD and hadronic sectors. 
The relation between a T and B obtained by elimination of the hadron size parameter agrees 
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very well with experiment. Starting from two experimental energies as input (541 and 1800 
GeV), this expression allows a determination of the correlation length and gluon condensate 
from high-energy data alone and gives good prediction of the remaining data. The results 
obtained are in good agreement with the correlations between the two QCD parameters 
obtained independently in lattice calculations and in the application of the stochastic vacuum 
model to hadronic spectroscopy. 

The present calculation is restricted to the lowest order non- vanishing contribution in the 
expansion of the exponential with functional integrations, and we may conclude from our 
results that this is justified for the evaluations of total cross-section and slope parameter. 
The resulting amplitude is purely imaginary and the p-parameter (the ratio of the real to 
the imaginary parts of the elastic scattering amplitude) can only be described if we go one 
further order in the contributions to the correlator. 

It is interesting to compare our results with Regge phenomenology, where the relation 
between the observables is given by cr^ egge = e °- 1616 ( B - B o) ^ obtained from a Regge am- 
plitude using the slope of the pomeron trajectory a / (0) pom = 0.25 GeV -2 . This relation 
requires an input pair , B at a chosen energy. Using as input the y/s = 541 GeV data 
Oq = 62.20 mb, and B = 15.52 GeV" 2 the line passes close to the CDF experimental 
point, instead of the E-710 point. If instead one uses as input Oq and B from the E-710 
experiment at y/s = 1800 GeV, the line shows a marked deviation at 541 GeV. This is rather 
intriguing, as it implies that the Regge formula favors the CDF experimental results at 1800 
GeV. 
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CAPTIONS FOR THE FIGURES 

Fig.l - Geometrical variables of the transverse plane, which enter in the calculation of the 
eikonal function for meson-meson scattering. The points C\ and C2 are the meson centres. 
In the integration, P2 runs along the vector Q(2, 1), changing the length z, which is the 
argument of the characteristic correlator function. In analogous terms, points Pi, Pi and 
Pi run along Q(l, 1), Q(l,2) and Q(2,2). This explains the four terms that appear inside 
the brackets multiplying k in the expression for the loop-loop amplitude. The length z' of 
the dot-dashed line is the argument of the Bessel function arising from the non-confining 
correlator Di; there are four such terms, appearing inside the brackets multiplying (1 — k) 



Fig. 2 - Relation between the two experimental quantities of the pp and pp systems. 
The values of a T at energies up to 62.3 GeV shown in this figure are the crJ om values as 
given by the parametrization aJ om = (21.70 mb)s a0808 of Donnachie-Landshoff. The values 
taken for B at the ISR energies are those for the pp system. The solid line represents 
B = B>a + Ca(ct t ) a , with values for A, and Ca obtained by fitting the data. The 
point at 1800 GeV used in the determination of the line is taken from the Fermilab E-710 
experiment, but the CDF point is also shown. 

Fig. 3 - Constraints on the values of (g 2 FF) and of the correlation length a. The solid 
line is the fit of our correlator to the lattice calculation^ as given in eq. flTB|) . The dashed 
line plots eq.(0), with p = 0.16 GeV 2 . The cross centered at a = 0.32 fm, < g 2 FF >=2.7 
GeV 4 shows our results. 

Fig. 4 - Energy dependence of the proton radius. The marked points are obtained from 
the total cross-section data (at the ISR energies the total cross-sections are represented by 
the pomeron exchange contributions). The two representations for the radius dependence 
are indistinguishable with the present data, but give quite different predictions for the cross- 
section values at the LHC energies. 
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TABLE AND TABLE CAPTION 
Table 1 . Values of the parameters for eq. ([12]) for k = 3/4 and k = 



1 



K 


a x 10 2 


13 


V 


7 


5 


6/13 


3/4 


0.6532 


2.791 


2.030 


0.3293 


2.126 


0.762 


1.0 


0.6717 


3.029 


1.859 


0.3118 


2.183 


0.721 
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proton e.m. radius = 0.862 fm 



solid line: S p =0. 572 + 0. 123 (logVs) 72 (fm) 
dashed line: S=0. 67 1 +0.057 logV^ (fm) 
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